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I. DEFINITIONS. 

The purpose, in the present discussion, is investigat- 
ing the possibility of an extension of Thermodynamics to 
systems which are in a stationary state but are subject to 
the action of conservative and nonconservative positional 
forces f and (therefore) also to the action of the forces 

i-pos v ' 

$ necessary to take away the heat thus generated. 

We first consider systems for which a finite microscopic 
mechanical model exists and are, therefore, described by 
equations of the form 

m'i = f_ po8 (x) + ±{x, x) d = F(x, x) 

and x is & point in an appropriate finite dimensional 
phase space (typically of very large dimension). If the 
force f is conservative then no thermostat is needed 

— pos 

and we suppose $ = (for simplicity). In general we 
call the force law $ a mechanical thermostat. 

(•) A key notion will be the phase space contraction rate 
a(x, x) which is defined as minus the divergence of the 
equations of motion: 

o~{x, x) = - d ±a F a (x, x) 

a=l 

(•) An equilibrium state will be a stationary probability 
distribution given by a density (one says a stationary 
11 absolutely continuous" distribution) on the phase space 
of a system which is subject only to conservative forces. 
We also identify the distribution with any point which is 
typical with respect to it: by typical we mean that the 
time averages of observables evaluated on the trajectory 
of the point are the same as the averages with respect to 
the distribution. 

(•) We suppose (ergodicity) that the time averages of ob- 
servables (just of the few physically relevant for macro- 
scopic Physics) are computable from any of the (equiva- 
lent) statistical ensembles: like the microcanonical en- 
semble. Hence an equilibrium state is identified with 
a probability distribution on phase space. A "typical" 
microscopic configuration, i.e. an initial datum in phase 
space which is not in a set of "unlucky cases" which, how- 
ever, form a zero volume set and are therefore (believed 



to be) 1 unobservable, will evolve in time so that the time 
averages of the observables (at least the few relevant for 
macroscopic Physics) are computable by the equilibrium 
state which has the correct values of the macroscopic pa- 
rameters: e.g. in the microcanonical case the energy U 
and the container volume V. When we speak of prop- 
erties of a single (typical) point in phase space, like of 
its "entropy" , we always mean the same property of the 
equilibrium state for which the datum is a typical one. 

Remarks: This already might be controversial: in fact the 
above (admittedly unconventional) definition of entropy 
has the following implications. A rarefied gas which ini- 
tially happens to have all molecules located in the left 
half of a container, because a separation wall has just 
been removed, setting the gas in macroscopic motion and 
out of the previous equilibrium state, will be an initial 
datum in phase space whose entropy is that of the same 
gas occupying the entire container and at the same tem- 
perature 2 (the difference residing only in the different 
dynamics that follows the removal the wall in the mid- 
dle of the container). Since a physicist would apply the 
Boltzmann equation to describe the evolution, the ques- 
tion arises about which is the place here of Boltzmann's 
H function, which is different if evaluated for the initial 
datum or for a datum into which the initial one evolves 
after a moderately large time (and in both cases it equals 
the classical thermodynamic entropy of the initial and fi- 
nal equilibria). 

(•) In trying to study nonequilibrium cases a conceptual 
difficulty must be met: if a system is subject to external 
non conservative forces then the thermostatting forces 
will have a non zero divergence and volume in phase space 
will not be preserved. The key idea (due to Ruelle)[l] is 
that in this case (as well as in the previous case) the 
nonequilibrium states will be the stationary states which 
are generated by time averaging of initial states that have 
a density or more satisfactorily, perhaps, by time aver- 
aging on the evolution of initial data that are typical 



1 This is not to be taken for granted, even though it is very often 
considered so: here I do not enter into discussing this mysterious 
assumption (as I have nothing to say). 

2 Assuming the gas to be ideal: we think here of Joule's experi- 
ment. In the following the temperature will have the dimension 
of energy, i.e. we call temperature T what is usually called kgT 
with kg being Boltzmann's constant 
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for probability distributions given by some (arbitrary) 
density. The latter states are called SRB distributions 
( Sinai, Ruelle,Bowen) [1-4]. Their appearance is natural 
since there will be no state (i.e. no probability distri- 
bution) which is stationary and at the same time is also 
given by a density. As in the equilibrium state we shall at- 
tribute to each individual point in phase space the macro- 
scopic properties of the stationary state which allows us 
to compute the averages of macroscopic observables on 
the motion of the given point. 

Remark: Systems that are chaotic in a mathematical 
sense (hyperbolic systems) can be shown, on rather gen- 
eral grounds, to have the property that there is only one 
SRB distribution (with the correct values of the macro- 
scopic parameters). [1,2,4] Therefore such systems, the 
only ones for which the above definition has a strict 
mathematical sharpness, verify an extension of the er- 
godic hypothesis: adopting the latter definition means 
believing that the system is chaotic enough so that typ- 
ical initial data generate a unique stationary distribu- 
tion with several features of the SRB distributions for 
hyperbolic systems: the assumption has been called 
chaotic hypothesis[5] and it represents (in our analysis) 
the noncquilibrium analogue of the classical ergodic hy- 
pothesis. 

(•) Given an initial state (a typical point or a distri- 
bution on phase space) it might be possible to define a 
function of it that will monotonically increase to a limit 
value which is the same for almost all data sampled with 
a distribution with a density on phase space (i.e. ab- 
solutely continuous). Since Boltzmann's H function is 
an example of such a function we shall call i7-functions 
all such functions. We recall that Boltzmann's H func- 
tion is defined by a coarse graining of phase space into 
"macrostates" determined by the occupation numbers 
f(p, q)d 3 pd 3 q of phase space cells d 3 pd 3 q around p, q and 
by defining H = — j f(p, q) log f(p, q) d 3 pd 3 q. Clearly 
here the cells size affects by an additive constant the 
actual value of H, which therefore should have no sig- 
nificance (at least in the classical mechanics context in 
which we are working) and only the variations of H can 
be meaningful. In a rarefied gas the Boltzmann equa- 
tion is believed to apply: [3] so that in the latter cases 
the function just defined is a nontrivial example of an 
ii-function. 

(•) Here we propose that even in the case of rarefied gases 
it is neither necessary nor useful that the H function is 
identified with the entropy: we want to consider it as a 
Lyapunov function whose role is to indicate which will be 
the final equilibrium state of an initial datum in phase 
space. This does not change, nor it affects, the impor- 
tance of Boltzmann's discovery that, in rarefied gases, 
the H function can be identified with the physical en- 
tropy whenever the latter is defined (i.e. in equilibria). 
Suppose that the initial state is chosen randomly with 
respect to a Gibbs distribution which is not the one that 



pertains to the given parameters that describe the system 
(e.g. volume V and energy U) but to other values: for 
instance it is chosen randomly with a Gibbs distribution 
t L u,v/2 that has the same energy but occupies half of the 
volume as in Joule's experiment. Then the inital and fi- 
nal values of the H function happen to coincide with the 
physical entropies of the Gibbs states nu,v and Hu,v/2 
(at least in a rarefied gas) which are given by the Gibbs' 
entropy. 3 

(•) The latter property could possibly be used to attempt 
a definition of entropy for states which are neither equi- 
librium nor stationary states: [6] however such a definition 
would be useful only in the special situations in which an 
H-theorem could be proved. That seems effectively to 
reduce the cases in which the notion would be useful to 
the ones in which an initial equilibrium state identified 
by some parameters (like U, V) evolves towards a final 
one identified by other values of the parameters. And 
even in such cases it is severely restricted to the rarefied 
gases evolutions in which the H-theorem can be proved: 
a statement that a model indpendent, universal, exten- 
sion of the above iJ-function, would have to be proved 
to have a monotonicity property, at least at a heuristic 
level, to avoid that its monotonicity becomes an a priori 
law of nature. 

(•) In general I would think that there will always be a 
Lyapunov function which describes the evolution of an 
initial state and is maximal on the stationary state that 
its evolution will eventually reach: however such a Lya- 
punov function may not have a universal form (unlike 
the H function in the rarefied gases cases) and it may 
depend on the particular way the system is driven by the 
external forces. After all the SRB distribution verifies 
a variational principle (Ruche) [1,2,4] which remarkably 
has the same form both in equilibrium and nonequilir- 
ium systems and one may imagine that in general it will 
be possible to define (on a case by case basis, I am afraid) 
a quantity that will tend to a maximum reached on the 
eventual stationary state. This picture seems to me sim- 
pler than trying to guess a (possibly nonexistent) general 



3 One should note that in the whole Boltzmann's work he has been 
really concerned with the approach to equilibrium: in our termi- 
nology he has been concerned with the problem of determining 
the stationary equilibrium distribution to which a given initial 
datum gives rise at large time: restricting consideration only to 
such cases one could well call the Boltzmann's //-function the 
"entropy" of the state as it evolves towards equilibrium. And, 
whatever name we give it, remains true that H is a measure of the 
disorder in the system. Our analysis here is intended to say that 
such an interpretation is not tenable when the system evolves 
towards a stationary noncquilibrium state. It must also to be 
said that even the H-theorem is not general because it applies 
only to rarefied gases, and even there it is an approximation: to 
generalize it to general situations, even when one only deals with 
approach to equilibrium, is a profound statement which should 
be substantiated by appropriate arguments. 
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definition of a quantity that would play the role of Boltz- 
mann's H. 



II. ENTROPY CREATION. 

The second law of Thermodynamics, in classical Ther- 
modynamic treatises, states: 

It is impossible to construct a device that, operating in 
a cycle, will produce no effect other than the transfer of 
heat from a cooler to a hotter body 

Of course this will be assumed to be a law of nature 
(Clausius) [7]. 

(•) The law implies that one can define an entropy func- 
tion S on all equilibrium states of a given system (charac- 
terized in simple bodies by energy U and available volume 
V) and if an equilibrium state 1 can be transformed into 
another equilibrium state 2 then 

using notations familiar from Thermodynamics: here the 
integral is over the transformation followed by the sys- 
tem in going from 1 to 2 and the /LQ is the heat that 
the system absorbs at temperature T from the outside 
reservoirs with which it happens to be in contact. The 
equality sign holds if the path followed is a reversible one. 

(•) One should note that the principle really says that 
the ^ does not depend on the path followed, if the 
path is a reversible sequence of equilibrium states; and its 
maximum value is reached along such paths. Existence of 
a path connecting 1 with 2 with S2 — Si < ^ would 
lead to a violation of the second law. 

(•) Is there an extension of the S2 — Si > j 1 rela- 
tion to noncquilibrium Thermodynamics? 4 Since there 
seems to be no agreement on the definition of entropy 
of a system which is in a stationary nonequilibrium and 
since there seems to be no necessity in Physics of such a 
notion, at least I see none, I shall not define entropy of 
stationary nonequilibrium systems (in fact the analysis 
that follows indicates that if one really insisted in defin- 
ing it then its natural value could, perhaps, be —00!). 

(•) In Thermodynamics one interprets — J x %jr as the 
entropy creation in the process leading from 1 to 2, be 
it reversible or not or through intermediate stationary 
states or not. The name is chosen because it is thought 
as the entropy increase of the heat reservoirs with which 



4 Which in a sense is tantamount of asking "is there a nonequilib- 
rium Thermodynamics?" 



the system is in contact and which are supposed to be 
systems in thermal equilibrium: so that their entropy 
variations are, in principle, well defined because they fall 
in the domain of equilibrium thermodynamics. In a gen- 
eral transformation from a state 1 to a state 2, both of 
which are stationary (non)cquilibrium states, following 
a path of (non)equilibrium stationary states (jP^ and in 
contact with purely mechanical thermostats one could 
consider the contribution to the entropy creation due to 
irreversibility in the process leading from 1 to 2 during a 
time interval [0, 0] to be 



A — c (&(x,x))^ t )dt = c / atdt 
Jo Jo 

where at d = (o~(x, is the average phase space con- 

traction computed in the state fi^\ This follows a re- 
cently proposed identification of a(x , x) as proportional 
to the entropy creation rate (here c is a proportional- 
ity constant) [8]. The quantity A is for mechanical ther- 
mostats the analogue of — j 1 %Q for the generic phc- 
nomenological thermostats characterized by a tempera- 
ture T. 

(•) In considering macroscopic systems one may imag- 
ine situations in which a system is partially thermostat- 
ted by mechanical forces for which a model considered 
physically reasonable is available 5 and partially by phe- 
nomenologically defined "heat reservoirs" characterized 
by a temperature T and able to cede to the system quan- 
tities fl Q of heat ( fl Q can have either sign or vanish) . 

In such more general settings a system in contact with 
several thermostats of which a few are modeled by me- 
chanical equations and a few others are unspecified and 
are just assumed to exchange quantities of heat jiQ at 
temperature T the second principle will be extended as 

assuming that the path leading from 1 to 2 consists en- 
tirely of nonequilibrium stationary states and that it lasts 
a time interval [0, 0]. Regarding the external thermostats 
as thermodynamic equilibrium systems — is the en- 
tropy increase of the reservoirs at temperature T and 
A = cf(a)t is interpreted as the entropy increase of the 
mechanical reservoir: if this interpretaton is accepted the 
above relation becomes the ordinary second law for the 
external reservoirs and could be read as "the entropy 



5 Typically these arc models of friction, as in the Navicr Stokes 
equation case in which the viscosity plays the role of a thermo- 
stat. Or in granular matter where the restitution coefficient in 
the collisions produces energy dissipation. Another well known 
example is in Drude's theory of electrical conduction. 
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of the rest of the universe does not decrease" (because 
A — j 1 %Q > o), where "universe" is not the astronomi- 
cal Universe but rather the collection of physical systems 
whose interaction with the system under study cannot 
be neglected. 

If a system is a in a stationary state in which at = 
(a) > (^-independent) this essentially forces us to say 
that its entropy, if one insisted in defining it at the time 



to of observation, could only be 
hinted above. 6 



- J_° oo (a) dt = — co as 



(•) Since the quantity a t is > (Ruelle)[2] under very 
general conditions (in fact always if the chaotic hypoth- 
esis is assumed to hold true) A > and the proposed 
estension is compatible with the main consequences of 
the second law. The constant c will be taken 1 because 
the factor 1 can be computed by studying the expression 
of a(x, x) in special models: at the moment, however, 
I see no immediate physical implications of the "univer- 
sality" of this choice of c and for the purposes of what 
follows c could be any constant, even non universal. 

(•) The above implies again that we shall not be able to 
define an entropy function unless a t = 0. The latter is 
the condition under which equilibrium Thermodynamics 
is set up: so that if one studies only transformations 
from equilibrium states to other equilibrium states it is 
possible to define not only the creation of entropy but 
the entropy itself (up to an additive constant). 

(•) A number of compatibility questions arise: suppose 
that the system evolves between 1 and 2 under the ac- 
tion of a thermostat which is modeled by forces that act 
on the system. For instance we can imagine a container 
with periodic bounday conditions, we call it a wire, con- 
taining a lattice of obstacles, which we call a crystal, and 
N particles, which we call electrons, interacting between 
each other and with the lattice via hard core interactions 
(say) and subject also to a constant force, which we call 
electromotive force, of intensity E; furthermore the par- 
ticles will be subject to a thermostat force of Gaussian 
type 7 (as it is essentially the case in Drude's theory as de- 
scribed in classical electromagnetism treatises [9]) which 
forces the particles to have an energy U — u(E) which 
is an assigned function of E. Suppose that the value of 
E changes in time (very slowly compared to the micro- 
scopic time scales) following a profile E(t) as drawn in 
Fig.l 



E(t) 



A 



Fig.l 
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In this case the wire performs a cycle which is irreversible 
and the integral j x 4jr- is because the system is adiabat- 
ically isolated (the thermostat being only of mechanical 
nature). The entropy variation of the system is defined 
because the initial and final state are equilibria and, since 
they are the same, it is 0: but there has been entropy cre- 
ation A > 0. 

(•) It is (perhaps) natural to define the "temperature" of 
a mechanical thermostat by remarking that in the mod- 
els studied in the literature it turns out that at is pro- 
portional to the work per unit time that the mechanical 
forces perform, the porportionality constant being in gen- 
eral a function of the point in phase space. Therefore we 
can call T ~ the time average of the proportionality con- 
stant between a and the work W per unit time that the 
mechanical thermostatting forces perform: in this way 
a t = 7^- and f® a t dt = J where d Q is the to- 
tal work 1 performed by the mechanical forces which we 
can (naturally) call the heat absorbed by the mechanical 



reservoir. 



(•) If one imagines that the above conducting wire model 
at the same time exchanges heat with two sources, ab- 
sorbing Q2 at temperature T 2 and ceding Qi at temper- 
ature Ti via some unspecified mechanism, and assuming 
that the profile of E{t) is as in Fig. 2 



Eq 



E(t) 



Fig.2 







where the value of E corresponds to a temperature T 
in the above sense, 
becomes 



The inequality -A + V*- < 



Jo 



0t T 2 T x T T 2 Ti ~ 



For instance, we see that if T\ = T 2 = T we have realized 
a cycle which is irreversible. In it a quantity of heat Q 2 — 
Qi — L, with L — Qo = TqA = Tq L a t dt is absorbed 
at a single temperature T = To and is transformed into 
the amount Q 2 — Q\—Loi work: however the inequality 
forbids this to be positive, as expected. 



6 If we imagine possible to replace a mechanical thermostat with 
a phcnomcnological thermostat at temperature T then the left 
hand side of the relation above remains unchanged but a part of 
—A becomes a contribution to the second addend. 

7 This is not the appropriate place to remind the Gauss' least 
constraint principle: it can be easily found in the literature [3]. 



III. MECHANICAL AND STOCHASTIC 
MODELS. 

A definition in Physics is interesting (only) if it is use- 
ful to decribe properties of the systems in which we are 
interested. Therefore having set the above definitions one 
should expect to be asked why all the work was made. 
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In this case the whole matter was originated by ef- 
forts to interpret results that started to appear in the 
late 1970's concerning numerical experiments in molecu- 
lar dynamics [8]. 

(•) It is obvious that in numerical experiments one needs 
to deal with a finite system (and even not too large): 
hence various models of thermostats were devised for the 
purpose of obtaining equations that could be transformed 
into numerical codes and studied on electronic machines. 
This was a theoretical innovation with respect to previous 
models which either relied on stochastic boundary con- 
ditions or, in the more sophisticated cases, with (poorly 
understood) systems with infinitely many particles. And 
it opened the way to import the knowledge in the theory 
of dynamical systems that had been being developed in 
the two preceding decades or so. 

The novelty with respect to stochastic thermostats was 
more conceptual than numerical. Given the number of 
particles a stochastic code is often only mildly more com- 
plex (and it could even be simpler) at least in the cases in 
which the noise is uncorrelated in time and acts on one 
particle at a time. This means that the resulting code 
does not require a longer running time than a determin- 
istic one: a fact that can also be seen by noting that a 
stochastic system can be regarded as a deterministic sys- 
tem with more degrees of freedom (i.e. the ones needed 
to describe the random numbers generators that one has 
to use and which, as it is well known, are simply suitably 
chaotic dynamical systems themselves). 

From the point of view of code writing this amounts 
at adding a few more particles to the system. 8 

(•) It is by no means clear that by using mechanical ther- 
mostats one can obtain physically realistic models nor, 
assuming that the stochastic dynamic models are more 
realistic, models behaving in as complex a way as the 
stochastic ones (typically consisting in boundary colli- 
sion laws in which the particles emerge with a maxwellian 
velocity distribution with suitable variance, i.e. suitable 
temperature). Understandably the matter is controver- 
sial but quite a few researchers think that this has been 
positively demonstrated by large amounts of work done 
in the last thirty years. 

Since stochastic models are just models with more de- 
grees of freedom it is tautological that there is equivalence 
between all possible stochastic models and all determinis- 
tic ones. The real question is whether the rather simple 
deterministic thermostats models that have been used 
are able to simulate accurately the stochastic models be- 



For instance if the stochastic thermostat is defined by requiring 
that upon collision with the boundary a particle rebounds with 
a maxwellian velocity distribution with dispersion (temperature) 
depending only on the boundary point hit then one needs three 
Gaussian random number generators, i.e. essentially three more 
degrees of freedom. 



lieved to be more realistic. 

(•) In my view it is likely that vast classes of thermostats, 
deterministic and stochastic as well, are equivalent in 
the sense that they produce motions which although 
very different when compared at equal initial conditions 
and at each time have, nevertheless, the same statistical 
properties [3,4]. And in my opinion there is already evi- 
dence that it is indeed possible to simulate the same sys- 
tem with simple deterministic thermostats or with some 
corresponding stochastic ones. 

Here the equivalence is intended in a sense that is fa- 
miliar in the theory of equilibrium ensembles: if one fixes 
suitably certain parameters then ensembles (i.e. time in- 
variant probability distributions in phase space) that are 
apparently very different (e.g. microcanonical and canon- 
ical) give, nevertheless, the same statistical properties to 
vast (not alt) classes of observables. If one fixes the en- 
ergy U and the volume V in a microcanonical ensemble 
or the inverse temperature and the volume V in the 
canonical ensemble then one obtains that local observ- 
ables have the same statistical distribution in the two 
cases provided the value of (3 is chosen such that the av- 
erage canonical energy is precisely U. 

One among the most striking examples of such equiva- 
lence (She, Jackson) [10] is the equivalence between the 
dissipative Navicr Stokes fluid and the Euler fluid in 
which the energy content of each shell of wave numbers 
is fixed (via Gauss' least constraint principle) to be equal 
to the value that Kolmogorov's theory predicts to be the 
energy content of each shell at a given (large) Reynolds 
number. Here one compares two very different mechani- 
cal thermostats. A more general view on the equivalence 
between different thermostats has been developed since. 
In fact many instances in which Physicists say that "an 
approximation is reasonable" really correspond to equiv- 
alence statements about certain properties of different 
theories (and in the best cases the statements can be 
translated into proper mathematical conjectures). 

(•) Coming therefore to consider more closely mechan- 
ical thermostat models the phase space contraction has 
turned out in many cases to be an interesting quantity 
often intepretable as the ratio between the work done 
by the thermostats on the systems and some kinetic en- 
ergy average: this led since the beginning to identify the 
phase space contraction rate with the entropy creation 
rate. The above "philosophical" considerations have been 
developed to give some background interpretation to the 
vast phenomenology generated by the new electronic ma- 
chines used as tools to investigate complex systems evo- 
lutions. 

A collision between the previously held views, egre- 
giously summarized by the book of De Groot-Mazur[ll], 
based on continuum mechanics and the new approaches 
based on transistors, chips and dynamical systems theory 
ensued: often showing that the two communities give the 
impression of not really meditating on each other argu- 
ments. 
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(•) Setting aside controversies it is interesting that the 
mechanical thermostats approach has nevertheless led to 
a new perspective and to a few new results. Here I men- 
tion the fluctuation relation: the phase space contrac- 
tion a(x(t),x(t)) which in various models has interesting 
physical meaning (like being related to conductivity or 
viscosity) is a fluctuating quantity as time goes on. Its 
average value in a time interval of size r divided by its 
infinite time average in the future (a) is a quantity p that 
still fluctuates. Of course it fluctuates less and less the 
larger is t and its probability distribution (easily ana- 
lyzable by observing it for a long time and by dividing 
the time into intervals of size r and forming a histogram 
of the values thus observed) is expected on rather gen- 
eral grounds to be proportional to e T ^ p ) for r large with 
£(p) being a function with a maximum at p = 1 (i.e. at 
the average, hence most probable, value of p) provided 
(a) > 0, i.e. provided there is dissipation in the system 
and the system is, therefore, out of equilibrium. 

If the dynamical equations are reversible, i.e. if there 
is an isometry I of phase space which anticommutes with 
the time evolution (x,x) — ► St(x,x) = (x(t),x(t)) in the 
sense that ISt — S-tl and furthermore I 2 = 1 then, 
provided the system motion is "very chaotic" , it follows 
that 

a-p) = a?) - p{<?} 

This is a parameter free symmetry relation that was 
discovered in a numerical experiment (Evans, Cohen, 
Morriss)[12] and which has been checked in many cases. 
By very chaotic one means that the motion of the system 
can be assimilated to that of a suitable Anosov flow whose 
trajectories fill densely phase space (transitive Anosov 
flow). 

(•) Indeed for transitive reversible Anosov systems the 
above relation holds theorem [5]. Since models of 
physical systems are not Anosov systems from a strict 
mathematical point of view the above relation cannot be 
applied, not even to cases in which the model is reversible 
and the trajectories are dense on the allowed phase space: 
the chaotic hypothesis says that the fact that the system 
is not mathematically an Anosov system is not relevant 
for physical observations, in most cases. This is similar 
to the statement that in equilibrium systems the lack 
of ergodicity of motions is irrelevant in most cases and 
averages can be computed by assuming ergodicity (i.e. 
by using the microcanonical distribution). 

If this is correct the above relation should hold: a non 
trivial fact to check due to the difficulty of observing such 
large fluctuations and to the lack of free parameters to 
fit the data, once they have been laboriously obtained. 

(•) When the forcing of the system is let to the above 



relation degenerates: not only (a) — > but also p itself 
becomes ill defined as its definition involves division by 
(a). Nevertheless by extracting the leading behavior of 
both sides the fluctuation relation leads to relations be- 
tween average values of derivatives of dynamical quanti- 
ties with respect to the intensity of the forcing, evaluated 
at zero forcing, and such relations can be intepreted as 
Onsager reciprocity relations and Green-Kubo expres- 
sions for suitably defined transport coefficients [4]. 

(•) Clearly a reversibility assumption on thermostats is 
a strong assumption and so is the chaotic hypothesis. 
Nevertheless the results are interesting and they seem to 
be among the few that can be obtained in a field which 
is well known for its imperviousness. The philosophical 
framework developed in Sections 1,2 helps keeping a uni- 
fied view on a subject that is being developed although, 
strictly speaking, one could dispense with the philosoph- 
ical view and concentrate on obtaining results that can 
be drily stated without appealing to entropy, entropy cre- 
ation, thermostats etc. 

(•) And one can go beyond various assumptions via the 
use of equivalence conjectures between different thcr- 
motats: for instance Drude's thermostat model which 
strictly speaking is not reversible is conjectured to 
be equivalent to a Gaussian thermostat which is re- 
versible. The Navier Stokes equation for incompressible 
fluids, clearly irreversible, is conjectured to be equiv- 
alent to a similar reversible equation [4] as the quoted 
experiment [10] shows and as other successive experi- 
ments seem to confirm[13]. The research along the just 
mentioned lines seems to go quite far and to lead not 
only to new perspectives but also to new results or con- 
firmations (i.e. non contradictions) of the general views 
in Sections 1,2. Doubts about the whole approach can be 
legitimately raised, and have been raised, on the grounds 
that the results are too few and too meager to be re- 
ally interesting: for instance one can hold against their 
consideration that they are not even sufficient to give 
some hint at a derivation of "elementary" relations like 
Fourier's law or Ohm's law. One can only say that time 
is not yet ripe to see whether the new methods and ideas 
lead really anywhere or at least to a better understand- 
ing of some of the problems that also the old ones have 
not been able to tackle, so far, (like the heat conduction 
laws or the electric conduction laws) in spite of intense 
research efforts. 

The above comments have been stimulated by continu- 
ous heated discussions held at Rutgers University in the 
course of the last few years: involving among many oth- 
ers S.Goldstein, J.Lebowitz, D.Ruelle. 

e-mail gallavotti@romal . inf n. it 

Piscataway, NJ, november 2002 
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